Quasiparticle cooling of a single-Cooper-pair-transistor 
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A superconducting tunnel junction is used to directly extract quasiparticles from one of the leads 
of a single-Cooper-pair-transistor. The consequent reduction in quasiparticle density causes a lower 
rate of quasiparticle tunneling onto the device. This rate is directly measured by radio- frequency 
reflectometry. Local cooling may be of direct benefit in reducing the effect of quasiparticles on 
coherent superconducting nanostructures. 



Coherent quantum nanostructures are highly sensitive 
to their thermal environment. In particular, supercon- 
ducting single charge devices are strongly affected by heat 
in the form of quasiparticles. The quasiparticle has a 
pronounced effect on these devices due to its electronic 
charge [1]. If a quasiparticle tunnels from the leads onto 
the device island the electrostatic energy of the system 
changes, and Cooper pair coherence is destroyed [I, 
Due to the detrimental effect on Cooper-pair coherence 
this effect is often known as 'quasiparticle poisoning'. Re- 
cent experiments have used high-bandwidth techniques 
to determine quasiparticle tunneling rates [4, 5, 6]. 

The temperature of a nanostructure can be re- 
duced by using on-chip, electronic refrigeration 
Superconductor-insulator-normal (SIN) tunnel junctions 
have been widely used for this purpose. This technique 
has been demonstrated to cool both metal islands [1, Q 
and suspended dielectric membranes |Tl|. In this 
Letter a reduction in quasiparticle poisoning is demon- 
strated by extracting quasiparticles from the supercon- 
ducting leads of a single-Cooper-pair-transistor (SCPT). 

The nanostructure consists of a single-Cooper-pair- 
transistor (SCPT) and a S1IS2 junction (Fig. 1(a)). 
The SCPT, using the quasiparticle poisoning effect, en- 
ables relative measurements of quasiparticle density in 
its leads. The S1IS2 junction acts as a way to reduce 
the quasiparticle population. Quasiparticle density de- 
creases exponentially with superconducting gap {Nqp ex 



vTAexp— ■^). Therefore, when the singularities in the 
density-of-states are aligned by applying a bias iVs^is^) 
to the S1IS2 junction such that {eVs-^is^ — ^l^i ~ ^21), 
the tunnel rate from the lower gap to the higher gap ma- 
terial is greater than in the reverse direction (Fig 1(c)) 
[l2:j . This allows cooling of the lower-gap region and the 
use of a S1IS2 junction as a refrigerator [l^. 

The device is fabricated by double angle evaporation 
of aluminium through a bilayer polymer resist mask. Be- 
tween the evaporation stages a controlled oxidation is 
performed to define the tunnel barriers. Thin aluminium 
films Jii which superconducting gap decreases with thick- 
ness , are used to generate the different gaps required 
for the Si IS 2 junction. To achieve continuous thin films 
the sample is placed on a low temperature stage (T ~ 200 
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FIG. 1: (Color online) (a) Schematic of the SCPT and S1IS2 
junction. The regions of two different gaps are indicated. 
Note that the lower gap region is isolated from the normal 
metal bond-pads (by the higher gap region) to avoid quasi- 
particle trapping, (b) Energy diagram of the SCPT, showing 
the potential minima for quasiparticles on the device island, 
(c) When the S1IS2 junction is biased to eVsiiS2 — |Ai — A2I 
quasiparticle extraction occurs from the SCPT reservoir. 



K) during evaporation. The A 1^2 regions have thick- 
nesses of 30 nm and 10 nm respectively. From previ- 
ous measurements on SIS junctions it was found that 
Al = 209 ± 11 /ieV and A2 = 250 ± 15 /ieV [15]. Device 
resistances, measured at 4.2 K, were Rscpt = 17.7 
and RS1IS2 = 3.4 kQ. The critical current of the S1IS2 
junction was suppressed from a maximum of Ic = 69 nA 
to /c < 2 nA by using a SQUID geometry. 

Radio-frequency reflectometry is used to measure the 
SCPT [16]. The device is embedded in a LC circuit and 
placed at milliKelvin temperature in a dilution refriger- 
ator. The circuit is resonant at 310 MHz and consists 
of a 470 nH chip inductor and a parasitic capacitance 
to ground. The complex reflection coefficient of a probe 
carrier signal at resonance depends on the presence of 
a quasiparticle on the SCPT. The carrier signal voltage 
biases the SCPT and its power is set to -94 dBm, cho- 
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sen to maximize the signal to noise ratio. The reflected 
carrier signal is demodulated and time records are taken 
using an oscilloscope. Details of the apparatus have been 
previously published [5[. 










1.0 








"5 


0.5 


O 




LL 

a: 


0.0 



fa 















5 10 15 
Time (^s) 



2 4 
Time (|lis) 



FIG. 2: (Color online) (a) Unaveraged gate sweep for the 
SCPT. (b) A single time trace at odd-integer charge. The 
points are separated by a time interval of 10 ns. (c) Histogram 
of the time-trace - solid lines are fits to a Gaussian distribu- 
tion, (d) Example histogram showing distribution of even 
times (no quasiparticle on island). The solid line indicates a 
fit to an exponential giving a time constant of te = 1.7 /xs. (e) 
Example of an odd-time (quasiparticle on island) histogram 
showing dual-Poissonian distribution. The time constants are 
toi = 22 ns and to2 = 1.3 /is. 

In an unaveraged gate-sweep on the SCPT both the 
intrinsic behavior of the SCPT and the effect of quasipar- 
ticle tunneling are observed (Fig. 2(a)). In the absence 
of quasiparticles only the 2e periodic supercurrent oscil- 
lations would be present. With quasiparticles present in 
the leads, two level switching occurs on the supercur- 
rent peaks as quasiparticles tunnel on and then off the 
island (Fig. 1(b)). A single shot time trace at the super- 
current maximum shows the bandwidth available with 
radio- frequency reflectometry (Fig. 2(b)). A rise time 
of approximately 30 ns is seen for a quasiparticle tunnel- 
ing event. A histogram of the time trace gives a dual 
peak distribution (Fig. 2(c)) conflrming that it is a two- 
level system. A real-time charge sensitivity is deduced, 
using dQ = ^j^-^^ ^ [17]. From the analysis described 
below the effective bandwidth of the measurement setup 



is determined to be 5 MHz. The charge sensitivity is 
dQ = 3.9 X 10~^ eHz~^-^, a value similar to previous 
frequency domain measurements. 

The time-traces are converted into a digital signal by 
comparing the recorded data with their median value. 
The times spent in the even (no quasiparticle) and odd 
(quasiparticle) states are then extracted and plotted in 
histograms. For the even times, the histogram is well- 
fltted by a single exponential (Fig. 2(d)) with time con- 
stant te- However, for the odd times, there are two time 
constants (Fig. 2(e)). This bi-exponential distribution is 
expected, due to the existence of both elastic and inelas- 
tic tunneling processes, and has been observed in recent 
measurements on Cooper pair boxes [6|. The effect of 
flnite measurement bandwidth is also apparent. It ap- 
pears as the peak at short times in the histogram. The 
maxima occurs at 200 ns (Fig. 2(e)), indicating a system 
bandwidth of 5 MHz [Ti. 
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FIG. 3: (Color online) (a) The extracted even time constant 
te, normalised to te{Vds = 0) = 1.18 /is, as a function of 
S1IS2 junction bias. There is a factor of two increase when 
eVsiiS2 = Ai — A2 indicating cooling of one of the reser- 
voirs. A simulation, based on balancing cooling power and 
heat transfer from the phonon system, is also shown, (b) 
The behavior of the odd time constant to2, normalised to 
to2{Vds = 0) = 58 ns under the same bias conditions. 

The time constant tg is proportional to the quasiparti- 
cle density in the leads Q . This allows relative changes in 
quasiparticle density to be determined. Due to the nature 
of the technique there is no, in principle, lower bound on 
the quasiparticle density that can be measured. At an ex- 
perimental temperature of 250 mK the quasiparticle den- 
sity in the leads is calculated to be 2.4 x 10^^ m~^, leading 
to an average number of quasiparticles in the Ai regions 
of each lead of ~ 1. The te time constant is measured as a 
function of S1IS2 junction bias (Fig. 3(a)). It is seen that 
te increases to a maximum at a voltage of Vs:^iS2 = 47 
/iV. The maximum, expected on the basis of the alu- 
minium fllm properties, eVs:^iS2 = A2 — Ai = 41 ± 18 
/ieV, shows reasonable agreement with the experimental 
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values. At the maximum tg is close to double its zero bias 
value indicating a two-fold reduction in quasiparticle den- 
sity. There is only a two-fold reduction since quasiparti- 
cles are only being extracted from one of the leads. From 
the closeness to a factor of two, the contribution from 
the drain lead to quasiparticle tunneling is inferred to 
become negligible. Either quasi-holes or quasi-electrons 
are extracted, depending on the polarity of the S1IS2 
junction bias. However the quasiparticle branches are 
strongly coupled, hence the effect on the time constants 
is symmetric with respect to bias polarity. The leads are 
separated from each other by the SCPT which regulates 
quasiparticle transfer by Coulomb blockade. As a result 
only a small effect of quasiparticle heat transfer through 
the SCPT is expected. Attempts to cool the source lead 
with an additional S1IS2 junction proved unsuccessful 
due to the radio- frequency signal being shorted out, this 
could be avoided by using a high value inductor (choke) 
on the source S1IS2 junction electrode to block the radio- 
frequency signal. This difficulty would also avoided in 
the case of a Cooper pair box which has only a single 
superconducting reservoir. The decrease in te to past it's 
zero-bias value, at Vds = 87 /iV, indicates an increase in 
quasiparticle density in the leads. This may be attributed 
to the onset of multi-particle tunneling processes [13]. 

A simulation of the change in te as a function of bias 
was performed. The cooling power is provided by the 
S1IS2 junction and is balanced by heat transfer from the 
phonon system which remains close to the lattice temper- 
ature. The expression for the heat transfer from phonons 
to quasiparticles is given by Qph-qp = T>V{T^^ - T^^) 
where S = 0.3 x 10^ Wm~^K~^ is the electron-phonon 
coupling constant for aluminium and 1/ = 2x0.1x0.03 
lirn? is the volume of the Ai region of the drain electrode 
[7|. In addition, heat transfer (Qs'i/S's) through a S1IS2 
junction can be numerically evaluated (Fig. 4(a)) [8I.[2q|. 
By consistently solving Qph-qp and Qsiis^ it is possible 
to determine a minimum quasiparticle temperature, and 
therefore a change in te- The result (Fig. 3(a)) signifi- 
cantly underestimates the observed effect. The assump- 
tion was made that the whole volume of the Ai region 
of the lead was cooled. However, one arm of the S1IS2 
junction is in proximity (~ 100 nm) to the SCPT poten- 
tially leading to a much stronger quasiparticle extraction 
effect near the measuring device (the SCPT). 

There is no significant change in to. As demonstrated 
in [13] , the effect of finite system bandwidth and the dis- 
crimination algorithm leads to an apparent change in one 
time constant as the other is varied. This can be cor- 
rected for but, since te is the more interesting time con- 
stant in this experiment, the un-processed values of the 
time constants are presented. 

Finally, the cooling power is plotted (Fig. 4(a)) as a 
function of bias across the S1IS2 junction, showing the 
divergence at eVs:^iS2 = — ^2!, when the singular- 

ities in the density of states are aligned. A simulation 
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FIG. 4: (Color online) (a) Simulated cooling power {QS1IS2) 
of a S1IS2 junction, at constant temperature, versus bias. 
Calculation parameters are similar to the experimental case: 
T=250 mK, Ai(2) = 210 (250) /xeV, R = 3.4 kQ. (b) The 
cooling power of the same junction biased close to — 
A2 — Ai as the temperature of the Ai material, corresponding 
to the lead of the SCPT, is reduced. 



was also performed of the cooling power, near it's max- 
imum value at Vs^iS2 — 40 /iV, as the temperature of 
Ai is reduced (Fig. 4(b)). The cooling power becomes 
negative for ^ = 0.86, close to ^ ~ showing the 
maximum temperature reduction achievable. For an es- 
timated initial temperature of 250 mK this would lead to 
a reduction in temperature to 215 mK and a consequent 
reduction in quasiparticle density, and increase in tg, by 
a factor of 5. To achieve a larger temperature reduction, 
A2 could be increased by using thinner aluminium films 
or other superconducting materials. Alternatively, alu- 
minium leads could be used with a lower gap material, 
such as titanium, for the island (therefore reducing Ai) 
[8]. 

In conclusion, this experiment demonstrates the po- 
tential of on-chip cooling to reduce the quasiparticle poi- 
soning effect in superconducting coherent devices. Fu- 
ture experiments on Cooper pair boxes should be able to 
achieve a dramatic reduction in quasiparticle poisoning 
by means of quasiparticle extraction. 
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